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A complete and exact solution of the ground-state problem for the Ising model on the Shastry- 
Sutherland lattice in the applied magnetic field is found. The magnetization plateau at the one third 
of the saturation value is shown to be the only possible fractional plateau in this model. However, 
stripe magnetic structures with magnetization 1/2 and 1/n (n > 3), observed in the rare-earth 
tetraborides RB4, occur at the boundaries of the three-dimensional regions of the ground-state phase 
diagram. These structures give rise to new magnetization plateaus if interactions of longer ranges 
are taken into account. For instance, an additional third-neighbor interaction is shown to produce 
a 1/2 plateau. The results obtained significantly refine the understanding of the magnetization 
process in RB4 compounds, especially in TmB4 and ErB4 which are strong Ising magnets. 



PACS numbers: 75.60.Ej, 05.50.+q, 75.10.Hk, 75.10.-b 

Geometric frustrations in lattice systems result in a 
rich variety of phenomena in both classical and quantum 
models. The investigation of such models, and even of 
their ground states is, however, a difficult problem. The 
first two-dimensional frustrated quantum model whose 
ground states have been determined exactly was intro- 
duced by Shastry and Sutherland in 1981 [l|. The 
Shastry-Sutherland (SS) lattice [Fig. 1(a)] is topolog- 
ically equivalent to the Archimedean lattice 3 2 .4.3.4 
[Fig. 1(b)]. In 1999, the SS model was shown to describe 
the magnetic properties of the compound SrCu2(BO,3) 
synthesized in 1991 Some time later, other quasi- 
two-dimensional compounds with magnetic atoms of each 
layer located on a lattice topologically equivalent to the 
SS one have been discovered. In particular this concerns 
the rare-earth tetraborides RB 4 (R = La-Lu) 0-0] ■ Some 
of these are regarded to be classical systems, since the 
magnetic moments of the magnetic ions are large. More- 
over, if the crystal field effects are strong enough, then 
the compounds can be described in terms of an effective 
spin-1/2 SS model under strong Ising anisotropy. This is 
the case of TmB 4 0, H| , ErB 4 0, Q , and H0B4 @ , where 
the easy-magnetization axis is normal to SS planes. 
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FIG. 1: (color online). (a) Shastry-Sutherland and (b) 
Archimedean 3 2 .4.3.4 lattices. They are topologically equiv- 
alent. 



SS magnets exhibit fascinating and puzzling sequences 
of fractional magnetization plateaus. For instance, 



plateaus at m/m s = 1/2,1/7,1/8,1/9,... of saturated 
magnetization m s have been observed in TmB 4 for tem- 
peratures below 4K when the field was normal to the 
SS planes; the 1/2 plateau is the major one. To explain 
the origin of this sequence of plateaus is a challenging 
task. Some efforts were made to do this through the 
Ising model on the SS lattice since TmB 4 is a strong Ising 
magnet as well as ErB 4 [5[ where a single 1/2 plateau 
was observed @[. Because of the strong frustrations, the 
solution of the ground-state problem for this model is 
difficult to obtain, therefore in Ref. [l| the Ising limit 
was analyzed only for the zero field. After thirty years of 
research, there is still no analytical solution for nonzero 
fields. There are only some numerical results showing the 
existence of a single fractional magnetization plateau at 
m/m s = 1/3. However, the results of numerical simu- 
lations cannot be considered as absolutely reliable since 
an inappropriate finite lattice sizes can lead to erroneous 
conclusions. For instance, in Ref. 0] a single magneti- 
zation plateau at m/m s = 1/2 was found for a system 
of 16 spins. More precise calculations, however, did not 
confirm this conclusion 

Here, we determine a complete and exact solution of 
the ground-state problem for the Ising model on the SS 
lattice u sing our own method that has been elaborated 
recently |lll [l2[ • We rigorously prove the existence of a 
single 1/3 plateau. As to magnetic structures with other 
fractional values of m/m s , they exist at the boundaries 
of full-dimensional regions of the ground-state phase di- 
agram. (A region in the parameter space of a model is 
full-dimensional if its dimension is equal to the dimension 
of the space). However, as the range of interactions in- 
creases, these structures become full-dimensional one by 
one and, hence, give rise to new fractional magnetization 
plateaus. For example, we show that an additional third- 
neighbor interaction produces a 1/2 plateau. A similar 
result was obtained numerically in Ref. [HI but we have 
found three different phases with m/m s = 1/2; one of 



2 



TABLE I: Basic rays and basic sets of the triangle configu- 
rations for the Ising model on the SS lattice in the magnetic 
field. 
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these is expected to be realized in ErB4. 

Let us consider the spin-1/2 Ising model on the SS 
lattice in the magnetic field 



H = J 1 £ 



ai °i ~ h J2^- (!) 



<y>i 



J 2 > rr ; n 

<ij>2 



Here Ci, <jj — ±1 are the spin variables; J\ and J2 are the 
interaction constants (for TH1B4 and ErB4, Jj « J 2 > 
0); < ij >i and < ij >2 denote the summation over 
edges and diagonals [not all, see Fig. 1(a)] of the squares, 
respectively; h is the magnetic field. 

The ground-state phase diagram for any Ising-type 
model is a set of convex polyhedral cones in the parame- 
ter space [Til] - (A polyhedral cone is the linear hull with 
nonnegative coefficients — the so-called conic hull — of a 
set of vectors. It is fully determined by its edges or vec- 
tors along them) . This follows from the obvious fact that 
the region corresponding to a ground-state structure is a 
solution of a set of linear homogeneous inequalities. The 
full-dimensional polyhedral cones fill the parameter space 
without gaps and overlaps. Herein we refer to a structure, 
which is a ground-state structure in a full-dimensional 
polyhedral cone, as "full-dimensional" and we refer to the 
correspon ding edg es (vectors) as the "basic ray" ("basic 



vectors") [ll|, [12J. The convexity property (the whole 
region in the parameter space, where a structure is the 
ground-state one, is always convex) makes it possible to 
find the ground-state structures in any point of the pa- 
rameter space if all basic rays are known as well as the 
ground-state structures in them. 

Let us consider the rays listed in Table I. We shall 
shortly see that these form a complete set of basic rays 
for the Ising model on the SS lattice. To determine the 
ground-state structures in these rays, let us rewrite the 
Hamiltonian (p} on the SS lattice as a single sum over all 



clusters in the form of a right-angled triangle with an SS 
diagonal being the hypotenuse, that is, 
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The numbers a^o, and a \% define a configuration 
(among the six unique possible ones) of the i-th trian- 
gle; <7io is the spin value in the vertex of the right angle. 
With each site being a vertex for three triangles, an ar- 
bitrary number a accounts for the fact that the energy 
contribution of a site can be distributed among these tri- 
angles in various ways. 

We say that a structure is generated by a set of triangle 
configurations if each triangle in the structure has a con- 
figuration belonging to the set. If, in a point (h, J\, J2), 
all configurations of the set have the same energy which 
is lower than the energies of all the remaining configura- 
tions, then the structures generated by the configurations 
of the set are ground-state ones in this point. Now, for 
each basic ray, we can find a set of triangle configura- 
tions which generate all ground-state structures in this 
ray. We call such sets "basic sets of configurations" [11 1 . 



They are given in the second column in Table I, and in 
the fourth column such a value of the "free" coefficient 
a is indicated, for which all configurations from the ba- 
sic set have the same energy Hi which is lower than the 
energies of all the remaining configurations. Herein solid 
and open circles denote spins up and down, respectively. 

Table I represents a complete solution of the ground- 
state problem for the Ising model on the SS lattice in 
the magnetic field. Using this table, we first of all have 
to determine the full-dimensional regions and structures. 
Here is an example sufficient to explain the way for do- 
ing this: The Neel structure is generated by the set of 
configurations 2b and oV This set is a subset of the ba- 
sic sets of configurations for ri, r 3 , r 5 , and rg" . Hence, 
the Neel structure is the ground-state one in these rays 
and, by virtue of the convexity property (that is, if a 
structure is a ground-state one in two points of the pa- 
rameter space, then this structure is a ground-state one 
in the entire line segment connecting these two points), 
in the whole polyhedral cone generated by these four vec- 
tors. A structure is full-, i.e., three-dimensional if it is 
generated by a set of triangle configurations that is a 
subset of at least three basic sets. All full-dimensional 
regions are listed in Table II. In its first, second, third, 
and fourth columns we see, respectively: (1) the notation 
of the region, (2) the set of triangle configurations gener- 
ating the ground-state structures in this region, (3) the 
basic rays which are edges of the region, and (4) the con- 
ditions for the existence of the region in the plane (h, J2)- 
The bar over the number of a region (structure) means 
that this region (structure) is symmetric to the region 
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FIG. 2: (color online). Stereoscopic three-dimensional 
ground-state phase diagram for the Ising model on the SS lat- 
tice (see also Tables I and II). Edges of the polyhedral cones 
(basic rays) are depicted in red; faces of cones are bordered by 
arcs of the unit sphere (heavy black lines). Unit circles (thin 
lines) in the three orthogonal coordinate planes are figured 
for better visualization. This figure is a 3-D image, therefore 
two parts of the figure are almost identical. 



TABLE II: Basic rays and generating triangle configura- 
tions for the full-dimensional ground-state regions of the Ising 
model on the SS lattice in the magnetic field. 
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(structure) with the same number but without bar with 
respect to the field inversion (the flip of all spins) . Now 
we can easily see that the vectors listed in Table I are 
basic vectors indeed and constitute a complete set. This 
is so, since the full-dimensional regions, determined from 
this set, fill the whole parameter space without gaps and 
overlaps (see the stereoscopic diagram in Fig. 2 and the 
two-dimensional diagram in Fig. 3). 

With the generating configurations for the full- 
dimensional ground-state structures being known, one 
can easily construct the latter. These are shown in Fig. 4, 
except for the fully polarized structure 1 (see also Fig. 3 
for the ground-state phase diagram). Phase 2 is a dis- 
ordered phase of Ising dimers (opposite spins on the SS 
diagonals). The disorder is two-dimensional, i.e., an ex- 
tensive entropy occurs. It can be easily shown than the 
entropy per site is equal to ks In 2/2 Phase 3 is the 
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FIG. 3: (color online). The two-dimensional ground-state 
phase diagram for the Ising model on the SS lattice (see also 
Tables I and II and Figs. 2 and 4). 
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FIG. 4: (color online). Full-dimensional ground-state struc- 
tures for the Ising model on the SS lattice in magnetic field. 
Unit cells of periodic structures are shown. Phases 2, 3, and 
4 are, respectively: the Ising dimer phase (opposite spins on 
the SS diagonals), the Neel phase, and the 1/3 plateau phase. 
Solid and open circles represent spins up and down, respec- 
tively. 



Neel phase. The magnetization of the phases 1, 2, 3, 
and 4 is equal to 1, 0, 0, and 1/3, respectively. Hence, 
there is a single fractional plateau at to/to, = 1/3, which 
confirms the numerical results of Refs. [a-llOl]. 

A very important advantage of our method is that it 
makes it possible to determine the ground-state struc- 
tures not only in full-dimensional regions but also at their 
boundaries. This makes it possible to analyze the ef- 
fects of longer-range interactions. The sets of generating 
configurations for the boundary structures are given in 
Table III. The sets are obtained as intersections of cor- 
responding basic sets. The structures observed in TmB 4 
emerge at the boundary between phases 3 (Neel phase) 
and 4 (1/3-plateau phase). They are generated by the 
configurations »b, 5b, and 0*. It is easy to verify that, in 
addition to structures 3 and 4, these configurations gen- 
erate a sequence of stripe structures as shown in Fig. 5 
and their mixtures. The stripes contain an even number 
2n of antiferromagnetic chains and are bordered by the 
ferromagnetic ones. We denote the structures composed 
of one type of stripe only by (3, 4)„ (n = 1 for structure 4 
and n = 00 for structure 3). The magnetization of struc- 
ture (3,4),i is equal to 2 n+i ■ ^ ne P er i°dic mixture with 
the smallest unit cell of consecutive structures (3, 4)„ and 
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TABLE III: Basic rays and triangle configurations for the two- 
dimensional faces of full-dimensional ground-state regions of 
the Ising model on the SS lattice in the magnetic field. 

Re- Basic Ground-state Existence in 
gions rays configurations the plane (h, J 2 ) 



1, 1 ri,r 2 

1, 2 r 2 ,r 4 

1, 3 n,r s 

1, 4 r 4 ,r 5 

2,4 r 3 ,r 4 

3, 4 r 3 ,r 5 



ob Sb 
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Ji < 0, J 2 < -2Ji 
Ji < 0, J 2 > -2Ji 
Ji > 0, J 2 < 
Ji > 0, J 2 > 
Ji > 0, J 2 > 2Ji 
Ji > 0, < J 2 < 2Ji 
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FIG. 5: (color online). Ground-state structures (3, 4) 2 and 
(3, 4)3 (m/m s — 1/5 and 1/7) at the boundary between phases 
3 and 4. The structures consist of ferro- (on yellow back- 
ground) and antiferromagnetic chains. Unit cells are indi- 
cated. 



(3,4) Il+ i (we denote it by (3,4) nj „ + i; see Fig. 6) has the 
magnetization equal to 2 „ 1 +2 • ^ n our future paper, we 
will rigorously show that structure (3,4)„ (n ^ 2) can 
become full-dimensional (i.e., it can produce a magneti- 
zation plateau) if the interaction range is not less than 
the distance between the successive ferromagnetic chains 
of the structure (3,4)„_i. At the boundary between re- 
gions 3 and 4, there can also occur multidomain struc- 
tures, similar to that shown in Fig. 7. They appear if 
some ferromagnetic chain forms a right angle. 

Let us consider some other boundaries. The struc- 
tures at the boundary between regions 1 and 3 (see Ta- 
ble III and Fig. 8) are constructed in the following way: 
Squares with SS diagonals should be chosen in such a 
manner that any couple of squares should have no com- 
mon sites. At these diagonals, the spins are orientated 
downwards and all the remaining spins are orientated up- 
wards. The structures at the boundary between regions 
1 and 4 (Fig. 8) are determined by a single condition: 
there cannot be two neighboring spins downwards nei- 



FIG. 6: (color online). Ground-state structure (3, 4) 2j 3 
(m/ms = 1/6) at the boundary between regions 3 and 4. 
This structure is the periodic mixture with the smallest unit 
cell of structure (3,4) 2 and (3,4) 3 (see Fig. 5). 



ther along SS diagonals nor along the edges. If an ad- 
ditional antiferromagnetic (as in TmB 4 ) third- neighbor 
interaction J3 along diagonals of "empty" (i.e., without 
SS diagonals) squares is included, then, from all possible 
structures at the boundary between regions 1 and 3 and 
1 and 4, only those structures become full-dimensional in 
which all the "empty" squares have the configuration *• 
(Fig. 8), since the energy should be as small as possible. 
These structures constitute disordered phases with mag- 
netization 1/2. One can show that the width of corre- 
sponding plateaus is equal to 4J3. Thus, contrary to the 
conclusion of Ref. [l3| , we can say that additional third- 
neighbor interaction is sufficient for the stabilization of a 
1/2 plateau. If J3 < 0, then among all structures at the 
boundary between phases 1 and 4 only the ordered struc- 
ture (1,4)_ shown in Fig. 8 becomes full-dimensional. It 
produces a 1/2 plateau as well. We suppose that this 
structure emerges in ErB 4 since J3 is expected to be neg- 
ative in this compound 5]. Interaction J3 lifts also the 
degeneracy of the Ising dimer phase. If J3 > (J3 < 0), 
then among all the structures of this phase, the structure 
with the minimum energy in the one for which spins at 
J3 bonds are different (identical). These structures are 
shown in Fig. 9. The Ising dimer phase exists just at the 
boundary between these. 
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FIG. 7: (color online). A multidomain structure at the bound- 
ary between regions 3 and 4 (m/m s = 0). Left and right 
ferromagnetic chains form right angles. 



Some ground-state structures, obtained here analyti- 
cally, coincide with those found numerically by different 
authors. Structure 4 was determined in Refs. [H [HI 13, 
[3| , structure 3 (Neel phase) in Refs. [li[l(i structures 
2 (Ising dimers) in Refs. [3, [1(3]. It seems that structures 
with magnetization 1/4, 1/5, and 1/6 coincide with those 
determined numerically in Ref. [14J for a spin-electron 
model (we cannot be sure, since the SS bonds are not 
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FIG. 8: (color online). Ground-state structures at the bound- 
aries between regions 1 and 3 and 1 and 4 (see also Table III). 
First and second figures present disordered structures. The 
right-hand parts of these figures show the arrangements of 
spins (still disordered) if an additional third-neighbor inter- 
action Jg > is included. All the "empty" squares have 
then the configuration *•• The third figure shows an ordered 
structure that becomes full-dimensional if J3 < 0. All the 
structures (1,3) + , (1,4) + , and (1,4)_ produce a 1/2 plateau. 
(See the text for details). 
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FIG. 9: (color online). Ground-state structures of the Ising 
model in the magnetic field on the SS lattice with additional 
third- neighbor interaction: (a) J3 > (collinear phase), (b) 
Js < (chessboard phase). The Ising dimer phase exists at 
the boundary between these two phases. 



indicated there). This is quite clear because the spin- 
electron interaction lifts the degeneracy at the boundary 
between regions 3 and 4 and some structures become full- 
dimensional. Our structures with magnetization » 1 . : 

n, 2n+l 

are different from those shown in Ref. [4{: all antifcr- 
romagnetic chains are shifted. The same concerns the 
structure with magnetization 1/3 shown in Ref. Q. 

The main conclusion of our study is that the Ising 
model with additional long-range interactions is sufficient 
to explain the origin of fractional magnetization plateaus 
in rear-earth tetraborides which are strong Ising mag- 
nets (TmB4 and ErB4). The long-range interactions are 
Ruderman-Kittel-Kasuya-Yosida ones |4l [l3j. since rear- 



earth tetraborides are good metals. Here we have studied 
the role of the first- and second-neighbor interactions and 
rigourously proved that they produce a single fractional 
plateau at 1/3 of saturated magnetization. We have also 
partially analyzed the effect of the third-neighbor inter- 
action and thus have shown that it produces at least five 
new phases; for three of these the magnetization is equal 
to 1/2. One of these 1/2-plateau phases is expected to 
emerge in ErB.4. Having analyzed the structures at the 
boundary between the Neel phase and the 1/3-plateau 
phase, we have also found the stripe structures with 
m/m s = 1/n (n 4); some of these have been observed 
in TD1B4. To finally explain the emergence of the se- 
quence of magnetization plateaus in T111B4, the effect of 
further-neighbor interactions should be investigated. 

We are grateful to T. Verkholyak for drawing our 
attention to the problem of fractional magnetization 
plateaus in SS compounds and for useful discussions. 
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